Entanglement Measure for Composite Systems 

V.I. Yukalov 
Bogolubov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research, Dubna 141980, Russia 



A general description of entanglement is suggested as an action realized by an arbitrary oper- 
ator over given disentangled states. The related entanglement measure is defined. Because of its 
generality, this definition can be employed for any physical systems, pure or mixed, equilibrium or 
nonequilibrium, and characterized by any type of operators, whether these are statistical operators, 

£f*l i field operators, spin operators, or anything else. Entanglement of any number of parts from their 

total ensemble forming a multiparticle composite system can be determined. Interplay between 
entanglement and ordering, occurring under phase transitions, is analysed by invoking the concept 

O^l , of operator order indices. 
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t-H ' Entanglement is the term Schrodinger coined for characterizing superposition of multipartite quantum states [1], 
which results in the appearance of specific quantum correlations between parts of a composite system. Nowadays there 
is a growing interest in studying entanglement due to its potential applications in quantum computing and quantum 
information processing [2,3]. In order to be a well defined characteristic, entanglement has to be quantifiable. One 
usually considers pairwise entanglement, for which several measures have been suggested, based on some kinds of 
reduced or relative entropies [2-8]. In the frame of a brief communication, it is impossible to present a detailed 
£C) • description of these measures, whose exhaustive account can be found in books [2,3], reviews [4-8], and references 
therein. But it is important to stress that the known entanglement measures are defined for quantifying only two- 
partite entanglement, and there is presently no definitive measure for entanglement between three or more subsystems. 
^ 1 Also, there is no a well defined entanglement measure for mixed multipartite systems. 

The aim of this communication is to introduce a general entanglement measure that would be valid for any system. 
Aiming at reaching a high level of generality, it is necessary, for a while, to leave aside all physical applications and 
C$ \ to focus our attention on the mathematical structure of the considered concept. 

Since entanglement deals with composite systems, we need, first of all, to concretize the meaning of the latter. 
In the present case, a system implies an object characterized by its space of states. Naming a space composite 
means that it is composed of some parts. Let the parts be labelled by an index i G I. The label manifold I can be 
discrete or continuous. In the simplest case, i = 1, 2, . . .. Let each part be characterized by a single-partite space Hi, 
which is a Hilbert space Hi = £{\ni >}, being a closed linear envelope of a single-partite basis {\m >}. Any vector 
ifii G Hi is presentable as an expansion ipi = J2 n a ni\ n i >■ The composite space H = ®{Hi is the tensor product, 
which is a closed linear envelope H — C{\{ni} >} of a multipartite basis >} whose vectors can be written as 

\{ni} >— (gii \m >■ For any ip G H, one has an expansion ip = J2{n } c {«;} \{ n i} >• 

Two remarks are in order. If the label manifold I is discrete, then ®{Hi is the standard tensor product with 
i = 1, 2, . . .. When I is continuous, then ®{Hi is the continuous tensor product, introduced by von Neumann [9] and 
employed for particular cases in Refs. [10,11]. Second, H is not compulsory a complete tensor product. It may happen 
that some selection rules are imposed on the latter, such as some symmetry requirements. Then H is a subspace of 
the complete tensor product and it is called [9] incomplete tensor product. 

Among all admissible vectors of H, one may separate out those of two types. One type forms the disentangled set 

T> = {®npi | ipi E Hi} (1) 

whose vectors / G T> C H have the structure of the tensor product / = ®j J2 n a m \ n i an( i are termed disentangled 
states. All other possible vectors of H constitute the complement H\T> whose elements cannot be presented as products 
of ipi G Hi and which are named entangled states. For instance, in the case of a bipartite system with two-dimensional 
single-partite spaces Hi, the examples of entangled states would be C12 1 12 > +C2i|21 > and cn|ll > +c 2 2|22 >, which, 
clearly, do not pertain to T>. 
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Entanglement, by its philological meaning, implies an action or a process by which disentangled states are trans- 
formed into entangled ones. A transforming action can always be described by an operator. Hence, one may consider 
entanglement produced by different operators. Let an operator A be given on TL. Acting on T>, this operator will, 
generally, transform disentangled into entangled states. We shall say that A is an entangling operator, provided that 
AT) = TL\T>. Of course, not each operator is entangling. And we shall call A a nonentangling operator when AV = T>. 

For each operator A onTL, we may put into correspondence a nonentangling operator A® having the structure of 
a tensor product <SnA\ of single-partite operators 

A\ = const Tr {Hjftt }A . (2) 

In order that the choice of a constant in Eq.(2) would preserve a scale- invariant form of A®, we require the validity 
of the normalization condition 



Tr n A = Tr v A® . (3) 
As a result, we obtain the nonentangling operator 



Trp ®i A\ 



A\ (4) 



associated with A. Note that the trace over the set T>, which is a restricted subset of the Hilbert space TL, is called 
the restricted trace. Such traces are widely used in statistical mechanics [10] and quantum information theory [2-7]. 
Their rigorous mathematical definition can be done by employing the corresponding projecting operators or, more 
generally, by invoking weighted Hilbert spaces [10]. In the present case, however, we do not need a general definition, 
since here the trace over V is applied only to the factor operators ®%A\, for which this reduces just to the shorthand 
notation Tr-p ®j A\ = J]; Tr Wi Ai. 

Now we come to the central problem of quantifying entanglement produced by an operator A on T>. The principal 
idea suggesting the way of constructing this measure stems from the following arguments. Entanglement, generally 
speaking, has to do with correlations between parts of a composite system. Interparticle correlations in physics are 
often connected with a kind of order classifying different thermodynamic phases and characterizing phase transitions. 
The level of ordering in physical systems can be described by order indices [12] advanced for reduced density matrices. 
This concept has been generalized by introducing operator order indices [13] defined for arbitrary operators. The 
definition of these order indices involves the norms of the corresponding operators. It is the operator norm that 
contains an essential information on the amount of order hidden in the action of this operator. Following this way of 
thinking, the amount of entanglement should also be related to operator norms. More precisely, we should correlate 
the actions on I? of a given operator A and of its nonentangling counterpart (4) by comparing the related norms 
\\A\\x> and ||A®||x>. Thus we finally arrive at the definition of the entanglement measure 

for the entanglement produced by an arbitrary operator A on the disentangled set V. As is evident, the measure 
e(A) = e(A,T>) is defined with respect to A as well as T>. But, for short, we may write e(A) when T> is fixed. The 
operator norms can be understood as those associated with the vector norms, so that \\A\\x> = sup^ eX) where 
||/||d = 1. The norm over a set V C TL is well defined [12], since it is straightforwardly reformulated to the norm 
||A||x> = WVvAVvWn over the Hilbert space TL by means of the projector Vt>, such that VvTL = T>. Though, in 
general, it is admissible to use different kinds of norms, everywhere in what follows the vector norms, associated with 
the related scalar products, are employed. This seems to be more convenient, in particular, because for maximally 
entangled two dimensional bipartite states, we get for the entanglement measure log 2. The base of logarithm may be 
any, though in information theory it is more customary to deal with logarithms to the base 2, when log 2 = 1. It is 
easy to show that the entanglement measure (5) possesses the following natural properties. 

1. Semipositivity: For any bounded operator A, 

s(A) > . (6) 

2. Nullification: Measure is zero for nonentangling operators having the structure of a tensor product A = A®, 

e(A 9 ) - . (7) 
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In particular, there is no sclf-cntanglement of a single part, when A = A\ = A®, and e{A\) = 0. The property (7) can 
be generalized to the case when A = @ v p v A® , where \\A® \\ v = \\A®\\ V and J2 V \P*\ = tnat is ' £ {®vPvA®) = 0. 

3. Additivity: For an operator A = ® V A„, 

e{® v A v ) =^e{A v ) . (8) 

V 

4. Invariance: Measure is invariant under local unitary operations 

e (® t U+A ® t Ui) = e(A) . (9) 



5. Continuity: If any considered operator A, being parameterized as A(t), with t G H., is continuous by norm, such 
that ||A(t)||x> — > ||A(0)||x) as t — > 0, then measure (5) is also continuous, 

e(A(t)) - e (A(0)) (t - 0) . (10) 

In this way, we may quantify entanglement produced by an arbitrary operator. Turning to physical systems, we 
could consider entanglement caused by any physical operator, for example, due to a Hamiltonian, number-of-particle 
operator, momentum, spin, and so on. It is, however, customary to examine entanglement only with respect to the 
von Neumann density operator p. The entanglement measure (5) can be easily calculated for this operator too, as we 
demonstrate below by several examples. For brevity, we shall write \i > instead of |rij >. 

(i) Einstein- Podolsky- Rosen states. The density operator for this famous example of a pure state is pepr = 
\EPR >< EPR\, where \EPR >= -^(\12 > ± |21 >). In order for readers to clearly understand how the measure 

is calculated, we illustrate for the present example all necessary details. Here the single-partite spaces are Hi, with 
i = 1,2. Each Hi is a two-dimensional Hilbert space, which is a span of the orthonormalized basis {\i >}. The 
composite space is H = Hi ®W 2 . The disentangled set V consists of the product functions / = ¥?i<S><^2, each ipi being 
a linear combination of the basis vectors \i >. The single-partite operator p\ = T^Hj^iPEPR, following definition 
(2), becomes p\ = |(|1 >< 1| + [2 >< 2|), given on Hi. The nonentangling operator (4) is p® — p\® p\. By their 
definition, all density operators are self-adjoint. For a self-adjoint operator p, the norm, associated with the related 
vector norm, can be written as ||p|| = su.pi i _f i i=i \{f,pf)\- In this way, ||pi||w 4 = sup J=12 < j\p\\j >= 1/2. Similarly, 
WpeprWv = sup itj < ij\pEPp\ij >= 1/2, while for the nonentangling operator, one has \\p®\\v = \\p\Wrd WpI \\n 2 = 
1/4. For measure (5), we immediately get s(pepr) — log 2, which is unity if the logarithm is to the base 2. 

(ii) Bell states. For the density operator ps = \B X B\, where \B >= ^ (|11 > ± |22 >), we again have 
e(p B ) = log 2. 

(iii) Greenberger-Horne-Zeilinger states. This is a three-partite state with pghz = \GHZ >< GHZ\, where 
\GHZ >= ^(|1H > ± |222 >). The corresponding measure is e{p GH z) = 21og2. 

(iv) Multicat states. These states are a generalization of the previous cases, the density operator being pmc — 
\MC X MC\, with \MC >= ci|ll . . . 1 > + c 2 |22 ...2 >, where |ci| 2 + |c 2 | 2 = 1 and N parts are assumed. The 
entanglement measure (5) is 

s(pmc) - (1- A01ogsup{|ci| 2 , |c 2 | 2 } . 
The maximum is reached for |ci| 2 = |c 2 | 2 = |, when s{pmc) — {N — 1) log 2. 

(v) Multimode states. Such states, that are a generalization of the multicat states, can be created in coherent systems 
of N parts, each of which can accept m different modes [14,15]. For the density operator pmm = \MM >< MM\, 
with | MM >= ^2 n c n \n . . . n >, where J2 n |c„| 2 = 1 and J2 n 1 = m, measure (5) becomes 

s(pmm) = (1 - N) log sup |c„| 2 . 

n 

Its maximal value happens for |c„| 2 = 1/m, when e(pmm) = {N — 1) logm. 

(vi) Hartree-Fock states. These describe N parts in different quantum states. The density operator is phf = 
\HF >< HF\, with \HF >= ^2 sym 1 12 . . . N >, where a symmetrized sum is implied, symmetric or antisymmet- 
ric, depending on the type of statistics, Bose or Fermi, respectively. Measure (5) takes the form 

N N 

e(pHF) = log — . 
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(vii) Mixed states. The von Neumann density operator can be of any type, but not only being related to pure states, 
as in the examples above. In general, it can be any nonequilibrium operator p(t), which leads to the evolutional 
entanglement, varying with time [15]. In the case of an equilibrium system, the von Neumann operator is p = 
Z~ 1 e~ f3H , where Z is a partition function, (3 is inverse temperature, and H is a Hamiltonian. Measure (5) provides 
the opportunity of quantifying entanglement for any p parts of an iV-partite system. To illustrate this, let us treat 
the case when a mixed p-partite state is obtained by tracing out N — p + 1 variables from the Hartree-Fock state of 
N parts, that is, p p HF = Tr« p+1 . . .Trn N Phf- The entanglement of this p-partite mixed state is measured as 

(N-p)\NP 
£(pH F )= l °g j^j • 

(viii) Statistical states. Mixed states in statistical mechanics, as is mentioned above, can be characterized by 
entanglement produced by the von Neumann statistical operator. However, this is not the sole possibility. Another 
way is to measure entanglement realized by reduced density matrices. This way in many cases may be much simpler, 
since the properties of the reduced density matrices have been thoroughly studied [12]. This method also provides 
a direct opportunity of quantifying entanglement for any p parts of a given statistical system. The procedure is as 
follows. Let x p = {x\, X2, ■ ■ ■ , x p } be a set of variables characterizing the annihilation, ip(x), and creation, ip^(x), field 
operators. A p-partite density matrix p p = {p p {x p ,x p )\ is defined as a matrix with respect to x p and xP , with the 
elements 

p p (x p ,x p ) = Tt^(x 1 )...^(x p )p^(x p )...^(x 1 ) , (11) 

where the trace is over the Fock space and p is a statistical operator. Under the variables x, one may mean, e.g., 
spatial coordinates, or momentum variables, or multi-indices in any convenient representation. The studied statistical 
system has N parts. These can be indistinguishable particles, with the field operators satisfying the boson or fermion 
commutation relations. Then a p-partite density matrix p p describes correlations between any p particles from the 
ensemble of N identical particles. It is reasonable to associate the single-partite space Hi with a span of the natural 
orbitals [12] that are the eigenvectors of the single-partite density matrix p\ = [pi(xi,Xi)]. In the considered case, for 
the nonentangling operator (4), we have 



=i P\ • (12) 



,p (N-p)\NP 

Keeping in mind identical particles, for which ||Pi||« 4 = ||jOi||wi, we obtain the entanglement measure (5) as 



(N-p)\NP\\p p \\ 



s{p p ) = log — , . (13) 



N\ \\ Pl 



What now is left is to find the norms of p p and p\ for a given statistical system, which can be done following the 
known prescriptions [12]. 

(ix) Spin states. Density matrices can be constructed not only of the field operators, as in Eq. (11), but of any 
other operators [13]. For instance, investigating spin systems, one may introduce spin density matrices [13] as follows. 
Let Si = {S" } be a spin operator associated with a lattice site i — 1,2, ... ,N. We may define a p-partite spin density 

matrix R p - ^ {ai/3} 



[p{a,/J 



as a matrix with respect to all indices, the matrix elements being 



R\$=TiSZ...SZpS%...Sf 1 1 , (14) 



with the trace over all spin states. The single-partite space Hi is defined as a span of the eigenvectors of the single- 
partite matrix i?i = [R°f]. The nonentangling matrix (4) is a p-fold tensor product R® = Ri ® R\ ® . . . ® The 
entanglement measure (5) becomes 

= 108 ■ (15 > 

(x) Phase transitions. Since entanglement is related to correlations existing in a composite system, it would not 
be surprising if entanglement would be sensitive to a physical order arising under phase transitions. Hence the latter 
can be accompanied by entanglement transitions. To prove this, let us consider some examples of phase transitions. 
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A. Bose-Einstein condensation. The changes in the reduced density matrices happening under this transition 
are well known [12,13]. Essentially above the condensation point, one has ||/9 p ||x> — HpiIIhi' because of which the 
entanglement measure (13) takes the form e(p p ) — \og(N — p)\N p /N\ typical of the mixed Hartree-Fock states. But 
below the condensation point, we have = N\/(N — p)\. Consequently, entanglement vanishes, e(p p ) = 0. 

B. Superconducting transition. Employing the properties of fermion density matrices, covered in great detail in 
book [12], we find the following. Above the critical point, entanglement measure e(p p ) is again of the Hartree-Fock 
form. But below the critical point, one has ||p P ||x> — c p Af( p_1 )/ 2 , when p is odd, and ||p P ||x> — c p N p / 2 , if p is even 
[12], where c p is a constant of order one. Thus, for measure (13), we obtain 



2 



logiV (p=l,3,...) 



| log TV (p = 2,4, 



where the large number of particles N ^> 1 is assumed. Consequently, entanglement increases, which is opposite to 
the case of Bose-Einstein condensation. 

C. Ferromagnetic transition. We shall study measure (15), based on spin density matrices [13]. For concreteness, 
let us keep in mind a Hcisenbcrg model with long-range interactions, when the mean-field treatment becomes asymp- 
totically exact in the thermodynamic limit. Then, using the properties of spin density matrices [13], we can derive 
measure (15). In the paramagnetic phase, we have 

but for the ferromagnetic phase we find e(R p ) = 0. Here the situation is analogous to Bose-Einstein condensation, 
where the arising long-range order leads to vanishing entanglement. Such a similarity can be understood if one 
remember that, under ferromagnetic phase transition, there occurs condensation of magnons [13]. 

As we see, phase transitions are really accompanied by entanglement transitions. However entanglement behaves 
differently under different phase transitions, sometimes vanishing but sometimes increasing, In order to fully under- 
stand the intimate relation between entanglement and ordering, occurring in physical systems, it is advantageous to 
resort to the notion of order indices that have been introduced for density matrices [12] and generalized to the case 
of arbitrary operators [13]. The operator order index for an operator A is 

lu(A) — l0 ^» 



log|Tr A\ ' 

For characterizing ordering in physical systems, the role of A is played by the appropriate density matrices, such as 
p p or R p . It is important that there may develop two types of long-range order, total and even [12,13]. When there 
is no any order, then \\A\\ <C |Tr A\ and u>(A) <C 1. If, under a phase transition, there develops total order, then 
\\A\\ ~ |Tr A\, so that the order index increases, lo(A) — > 1, but at the same time, because of the normalization 
condition (3), entanglement vanishes, e(A) — ► 0. This is the situation taking place at Bose-Einstein condensation or 
ferromagnetic transitions. Another case happens under the appearing even order, when ||p p || ~ ^/Tr p p /N, if p is 
odd, while \\p p \\ ~ WTrp p , if p is even. Then the order index increases to co(p p ) = (p — l)/2p, when p is odd, and 
to Lo(pp) = 1/2, if p is even. This results in the increase of the entanglement measure s(p p ), as it is shown above for 
superconducting transition. 

Concluding, a general entanglement measure is introduced, which describes entanglement realized by an arbitrary 
operator. For physical systems, it is convenient to consider entanglement caused by reduced density matrices. The 
relation between entanglement measure and order indices is investigated. 
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